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Abstract. We present mid IR spectro-photometric imag- 
ing of a sample of eight nearby (D < 240pc) Herbig 
Ae/Be stars. The spectra are dominated by photospheric 
emission (HR6000), featureless infrared excess emission 
(T Cha), broad silicate emission feature (HR5999) and the 
infrared emission bands (HD 97048, HD 97300, TY CrA, 
HD 176386). The spectrum of HD179218 shows both sili- 
cate emission and infrared emission bands (IEB). All stars 
of our sample where the spectrum is entirely dominated 
by IEB have an extended emission on scales of a few thou- 
sand AU (~ 10"). We verify the derived source extension 
found with ISOCAM by multi-aperture photometry with 
ISOPHT and compare our ISOCAM spectral photometry 
with ISOSWS spectra. 
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1. Introduction 

The mid- infrared spectrum of Herbig Ae/Be stars is rich in 
information that can be used to improve our understand- 
ing of the circumstellar environment. The wavelength in- 
terval between 6 and 13 fim contains most of the infrared 
emission bands (IEB), as well as the broad silicate feature 
centered at about 9.7 /xm. 

From ground-based and KAO spectroscopy (Roche et 
al. 1991, Schutte et al. 1990) it is known that IEB exist 
in the spectra of some Herbig Ae/Be stars. A large num- 
ber of Herbig Ae/Be stars have been measured with the 
Infrared Space Observatory (ISO). The ISO sample will 
allow us to determine the properties and nature of the 
IEB carriers with unprecedented accuracy (Waelkens et 
al. 1996). These data will be crucial in solving the puzzle 
of the mid-infrared energetics of Herbig Ae/Be stars, by 
determining to what extent IEB emission contributes to 
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the observed mid-infrared luminosity, and how much IR 
emission comes from a circumstellar disk in Herbig stars 
of different spectral type and age (Kenyon & Hartmann 
1991, Hillenbrand et al. 1992, Natta et al. 1993). 

A very important contribution to our understanding 
of the IEB and of their role in Herbig Ae/Be stars is pro- 
vided by spatial information on the extension and distri- 
bution of the emission at various wavelengths. For exam- 
ple, Prusti et al. (1994) were able to infer the presence of 
infrared emission bands (IEB) and to determine some of 
the properties of the distribution of the carriers around 
Herbig Ae/Be stars from multi-aperture photometry in 
the three narrow-band filters Nl, N2 and N3. Natta & 
Kriigel (1995), using the radiation transfer code developed 
by Siebenmorgen et al. (1992), computed the expected in- 
tensity profile in the IEB and in the adjacent continuum as 
a function of the spectral type of the exciting star and of 
the surrounding shell parameters. However, until recently 
the scarcity of observed intensity profiles has greatly lim- 
ited the usefulness of these calculations. 

Spectro-photometric imaging observations using the 
circular variable filters (CVF) of ISOCAM (Cesarsky et 
al. 1996) are capable of producing the spatial information 
we want. In addition to the determination of the source 
size at the various wavelengths, they can also reveal ad- 
ditional emission structures and deviations from spherical 
symmetry. We present in this paper the results of ISO- 
CAM spectro-photometric, ISOPHT (Lemke et al. 1996) 
multi-aperture and ISOSWS (de Graauw et al. 1996) grat- 
ing scan observations of a small sample of Herbig Ae/Be 
stars that we consider particularly suited to this purpose. 

Our sample includes 5 stars of spectral type B9/A0 
(HD 179218, HD 176386, TY CrA, HD 97048 and 
HD 97300). They cover a small age interval, from 1 Myr to 
> 3 Myr (Table 1). Some of these stars (i.e., TY CrA and 
HD 97048) are known to have strong IEB features (Roche 
et al. 1991, Deutsch et al. 1995). In addition, Prusti et al. 
(1994) could infer from multi-aperture study the presence 
of IEB for HD176393, HD97048 and HD 97300. 

In addition to these 5 stars, we have included in our 
sample two stars of later spectral type, namely HR 5999 
(spectral type A5/7) and T Cha (spectral type G8). The 
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aim was to check the dependence of IEB band shapes and 
intensity profiles on the spectral type of the exciting star. 
Ground-based observations in the 10 /im region (Wooden 
1994) seem to indicate that, as expected, IEB are weak or 
absent in Herbig Ae/Be stars of later spectral type (see 
also the 3.3 /im survey of Brooke et al. 1993). Much more 
surprisingly, it was also found that the silicate feature ap- 
pears in emission only in stars of later spectral types. 
There is no obvious reason why silicate emission should 
disappear in hotter stars. Its absence is probably due to 
the geometrical distribution of the dust around the star. 
If we believe that the silicate emission is linked to the ex- 
istence of a circumstellar disk (Chiang & Goldreich 1997), 
it is possible that early-type stars do not show the sili- 
cate feature because they do not have disks. This would 
confirm the results obtained with millimeter interferom- 
etry, namely that Herbig stars of spectral type A show 
evidence of disks while very few Herbig B stars do (Di 
Francesco et al. 1994). It is possible that silicate emission 
disappears in B stars as they age, dissipating their disks 
and pushing farther away the residual envelopes. The dis- 
covery of "transitional" objects, showing both IEB and 
silicate emission would be very interesting in this context. 

A convincing answer to these questions requires a sta- 
tistical (i.e. large) sample of spectra and intensity maps 
for Herbig Ae/Be stars of different spectral type and age. 
However, we feel that our results can provide interesting 
clues to a possible answer. 

This paper contains in Sect. 2 a description of the ob- 
servations and data reduction techniques. The results for 
the individual stars are presented and discussed in Sect. 3. 
A summary is given in Sect. 4. 

2. Observations 

Spectro-photometric images have been obtained with ISO- 
CAM for the stars in our sample in a number of CVF wave- 
length intervals, selected to probe the emission in IEB at 
6.2, "7.7", 8.6 and 11.3/^m. Additionally, one interval was 
chosen around the peak of the silicate feature at 9.7 /im. 
Typically, we have taken 25 exposures of 2.1 s at each CVF 
wavelength. We used a gain of 2 and the lens providing 3" 
pixel field of view. The basic data reduction steps are de- 
scribed in Siebenmorgen et al. (1999). After dark current 
subtraction, glitch removal and flux transient correction, 
the exposures are coadded in each CVF step. Photome- 
try is obtained applying the spectral response function. 
After background subtraction we simulate multi-aperture 
photometry on the images. 

In order to understand the photometric uncertainty of 
this procedure we analyzed data of the calibration star 
S Draconis and compared it with the Kurucz model. The 
largest uncertainty was found in the overlap region of 
CVF1 and CVF2 and is typically of the order of < 5%. The 
aperture sequence allows us to determine the full width at 
half maximum flux (FWHM). For a point source (6 Dra- 



conis, as shown in Fig. yj), we derive in the wavelength 
range considered, typical sizes of 3-4" (FWHM) with an 
uncertainty of about half a pixel size (~ 1.5"). 
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Fig. 1. The size (FWHM) of the comparison star S Dra 
as derived from multi-aperture photometry on the CVF 
images. 

We searched in the ISO archive and whenever possi- 
ble we compared our ISOCAM spectral photometry with 
ISOSWS data. The ISOSWS spectra are from the off-line 
processing version 8.6 (Leech et al. 2000). 

The CVF images are contaminated by ghosts and 
stray light. In order to independently check the FWHM 
as derived by ISOCAM we compare them to results ob- 
tained by ISOPHT. We selected multi-aperture photome- 
try (PHT04, Laureijs et al. 2000) in 11 apertures ranging 
from 5" up to 180" at an integration time of 32s per fil- 
ter. The dark subtracted, deglitched and linearised signals 
of our target stars are compared to an identical PHT04 
sequence measured on a point source (HR7127). Such a di- 
rect comparison is necessary because the point sources do 
not exactly follow the theoretical PSF behaviour (Miiller 
2000). However, adequate point source sequences exist 
only in a limited number of filters, in practice we can 
present our PHT04 measurements only at 7.3/^m. 

3. Results for Individual Stars 

The stars in our sample are listed in Table 1. Column 1 
gives the name of the star, Column 2 the spectral type, 
Column 3 the distance, Column 4 the star effective tem- 
perature, Column 5 the luminosity, Column 6 the age and 
in Column 7 we note the molecular cloud to which the 
star belongs. Distances, effective temperatures, luminosi- 
ties and ages are taken from van den Ancker et al. (1996). 
For TY CrA, for which no Hipparcos distance is avail- 
able, we have taken the same distance as of HD 176386 
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Name 


ST 


Distance 
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Age 


Cloud 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


HD 97048 


AO 


180l 3 2 ° 


10000 


41 


>2 


Chamaeleon I 


HD 97300 


B9 


188tf 


10715 


37 


>3 


Chamaeleon I 


HD 179218 


B9 


24411° 


10715 


316 


0.1 


no apparent cloud 


HD 176386 


B9IV 


1401^ 


10715 


49 


2 


Corona Australis 


TY CrA 


B9 


140 1 3 2 ° 


10715 


98 


1 


Corona Australis 


HR 5999 


A5-7III/IV 


210l 5 3 ° 


7943 


85 


0.5 


Lupus 3 


HR 6000 


B6 V 


241+11 


12882 


263 


>0.5 


Lupus 3 


T Cha 


G8 


66lg 


5888 


1.3 


>13 


Chamaeleon I/II 



and computed the luminosity from the observed V, (B-V) 
(Shevchenko et al. 1993), corrected for extinction with the 
assumption that the ratio of selective-to-total extinction 
is i?=5.1 (Bibo et al. 1992). The age of the star is then de- 
rived by comparing its location on the HR diagram to the 
evolutionary tracks of Palla & Stahler (1993). The same 
procedure has been used to estimate the luminosity and 
age of HR 6000. This procedure is similar to that used by 
van den Ancker et al. (1996). 

3.1. HD 97048 and HD 97300 

These two stars are located in the Chamaeleon I cloud, 
and have very similar properties. Their spectral type is AO 
(HD 97048) and B9 (HD 97300); they have luminosities of 
40 and 37 L Q , respectively and lie very close to the ZAMS. 
Both stars have been detected by IRAS. HD 97048 has a 
near-infrared excess, while HD 97300 has not. 

The mid-infrared spectrum of HD 97048 is known to 
be rich in IEB from ground-based and KAO observations 
(Roche et al. 1991, Schutte et al. 1990; Brooke et al. 
1993). Prusti et al. (1994) found that the emission was 
extended in both sources (by comparing aperture pho- 
tometry between 5" to 16"), but more so in HD 97300 
than in HD 97048. Our ISOCAM and a ISOSWS spec- 
trum of HD 97048 is shown in Fig. ^. It is indeed dom- 
inated by IEB emission centered at 6.2, "7.7", 8.6 and 
11.3/im. HD 97048 is extended on scales of a few arcsec- 
onds (5-10"). Consequently the ISOSWS spectrum (Van 
Kerckhoven et al. 1999), which is flux calibrated assum- 
ing a point source, underestimates the total emission; its 
high spectral resolution shows that the "7.7" /im band is 
due to two separate components, which are unresolved by 
ISOCAM. 

A deeper full CVF scan combined with photometric 
mid IR images of HD97300 has been presented by us in a 
separate paper (Siebenmorgen et al. 1998). It was found 
that the complete mid IR spectrum of the source is dom- 
inated by a huge ring structure of about 7500 AU in size. 
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Fig. 2. Top: Mid infrared spectrum of HD97048. The ISO- 
CAM CVF spectrum is indicated by O and connected by 
the solid line. The ISOSWS spectrum is shown by the 
dotted line. Bottom: The FWHM as derived from multi- 
aperture photometry on the CVF images. 

It was proven, by means of a dust model including tran- 
siently heated particles, that all the ring emission is due 
to the IEB carriers. The amount of circumstellar matter 
was found to be rather small (~ 0.03 M Q ). In Fig. || we 
compare the full CVF wavelength spectrum with the ob- 
serving strategy of sparse CVF scans, presented in this 
paper. One notices that the absolute photometry of both 
observations are identical down to a few %. Good agree- 
ment is also found for the derived FWHM of both CVF 
scans. We derive a FWHM of 13-19", which gives a clear 
sign of a huge extended nebulae as was independently ver- 
ified by stray light free, narrow band ISOCAM images. 

The FWHM of 14 ±3" as derived from our ISOPHT 
multi-aperture sequence at 7.3/im, presented in Fig. |J, is 
consistent with the CVF results. 
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Fig. 3. Top: Mid infrared spectrum of HD97300. The CVF 
full scan by Siebenmorgen et al. (1998) is shown as solid 
line and the CVF sparse scan of this paper by (O). Bot- 
tom: The FWHM as derived from multi-aperture photom- 
etry on the CVF full scan (line) and sparse scan (O). 
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Fig. 4. ISOPHT multi-aperture sequence on a point 
source (HR 7127) and the extended source HD 97300 at 
7.3 /im. The line indicate the background estimate and the 
intersection at 0" gives the source flux. The dashed line 
shows for HD 97300 the distance at half the source flux 
(FWHM). The FWHM of HD 97300 as determined from 
this figure is 14 ± 3". 



3.2. TY CrA and HD 176386 

The triple system TY CrA (spectral type B9 of the pri- 
mary) and the binary HD 176386 (spectral type B9 of the 
primary) are located in the Corona Australis molecular 
cloud. The stars are separated from each other by about 



55". Both stars are located very close to the ZAMS. Both 
have been detected by IRAS, but show no near-infrared 
excess, if anomalous extinction is taken into account (Bibo 
et al. 1992). TY CrA has been detected at 1.3 mm, but 
the amount of circumstellar matter associated to it is quite 
small (Hillenbrand et al. 1992; Natta et al. 1997). The mid- 
infrared spectrum of TY CrA was measured in a small 
aperture (4") by Roche et al. (1991), who found strong 
IEB. Multi-aperture photometry of HD 176386 by Prusti 
et al. (1994) indicate extended emission. 

In Fig. H we show the morphology at 11.28 /im. At 
about 8" SE from the peak we see a bar like structure 
(TY CrA bar). The bar is not caused by stray- light or 
ghost components, which can be predicted with the optical 
model of ISOCAM. The spectrum of the spherical compo- 
nent of TY CrA is shown in Fig. ^. Photometry was done 
in half circle apertures centered on the star at the oppo- 
site side of the bar. This flux was multiplied by two to get 
an estimate of the total flux without the bar component. 
This procedure gives an estimate of the flux in the spheri- 
cally symmetric envelope around the star. The spectrum is 
rich in IEB and has a very pronounced 8.6 /im band. The 
11.05 /im feature detected by Roche et al. (1991) is not re- 
solved by ISOCAM but clearly present in both ISOSWS 
spectra (Corporon et al. 1999). The source is very ex- 
tended (~ 2000 AU) in both the IEB and the continuum. 
Consequently the absolute flux in both ISOSWS spectra 
is underestimated by ~ 40%. 

In order to deduce the spectrum of the bar like com- 
ponent we have mirrored the pixels in the half circle used 
for 41" aperture photometry and subtracted these images 
from the original background cleaned maps. The spectrum 
of the TY CrA bar is shown in Fig. It shows again strong 
IEB but weaker " 7.7" and 8.6 /im bands as compared to the 
spherical component of TY CrA. 

As in TY CrA, we see two peaks of emission also in 
HD 176386. The maximum brightness is at the stellar po- 
sition. Photometry of the spherical component was done 
with half circles in a similar manner as for TY CrA. The 
spectrum is rich in IEB. The band ratios are quite similar 
to the bar component of TY CrA (Fig. |). The FWHM 
as derived from the ISOCAM images are indicative of a 
large extended halo. They are confirmed by the ISOPHOT 
multi-aperture sequence at 7.3/im. In Fig. 10 we show the 
ISOPHT measurements of HD176386. For apertures above 
80" the signals are strongly contaminated by the adjacent 
source TY CrA. We derive a source size of 9 ±3" (FWHM) 
for HD176386 and 13±3" (FWHM) for TY CrA (Fig. |ll|). 
Both size estimates are consistent with the CVF results. 

In order to derive photometry of the second component 
of HD176386, we mirrored the pixels of the half circle 41" 
aperture and subtracted the result maps from the original 
background cleaned images. The spectrum of the second 
peak, "HD176386 bar", is computed in a 14" aperture. 
The source is clearly detected at all wavelengths (Fig. [|). 
Its spectral shape, in particular of the 6.2/" 7.7" band ra- 
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tios, are quite different from the HD176386 main compo- 
nent (Fig. ^) but similar to TY CrA main component and 
HD97048 (Fig. |). 

3.3. HD 179218 

HD 179218 (also called MWC 614) is a young star of spec- 
tral type B9. It is listed as a Herbig Ae/Be star by The 
et al. (1994). It has no companion down to a distance of 
0.4" (Pirzkal et al. 1997). To the best of our knowledge, it 
has not been observed in the near-infrared, but has IRAS 
fluxes in the four bands and shows a 10 /im emission fea- 
ture in the low-resolution IRAS spectra. 

The multi-aperture CVF spectrum is shown in Fig. |l2| . 
There are pronounced IEB visible at 6.2, "7.7" and 8.6 /im 
on top of a strong silicate emission feature. As discussed 
by Waelkens et al. (1999) the silicate band of HD179218 
is dominated by substructures due to crystallinity. Those 
are absent in spectra of the interstellar medium or in the 
disks of the youngest stellar objects. We do not find any 
clear sign of an extended structure in the ISOCAM or 
in the ISOPHT observations. As expected for a point- 
like sources, the ISOSWS photometry is in line with our 
ISOCAM spectrum. 

3.1 HR 6000 and HR 5999 

HR6000 and HR5999 form the proper motion binary 
system Dunlop 199. They are located in the Lupus 3 
star forming region which is at a distance of 140±20pc 
(Hughes et al. 1993). Van den Ancker et al. (1996) classi- 
fied HR6000 as a young main sequence object with spec- 
tral type B6. The strong X-ray emission is accounted for 
by a T Tauri companion (Zinnecker & Preibisch 1994). We 
have a detection of HR 6000 above the background level. 
The star is close to the much brighter target HR 5999 so 
that only the central five pixels could be used to deduce 
the spectrum shown in Fig. [l3|. By comparing the spec- 
trum with a black body of 13000 K, we notice that it is 
featureless. 

The pre-main sequence star HR5999 (also called 
V856 Sco) has spectral type A7III (The et al. 1994) and 
has Rossiter 3930 as companion (Stecklum et al. 1995), 
which probably accounts for the strong X-ray emission. 
HR5999 has a strong near-IR excess (Hillenbrand et al. 
1992), and has been detected by IRAS at 12, 25 and 60 
/im, as well as at millimeter wavelengths (Henning et al. 
1994). The amount of circumstellar matter associated to 
the star is of the order of 0.006 M©. Our mid-infrared spec- 
trum of HR5999 (Fig. [l4]) shows clearly a broad silicate 
emission feature. In the CVF spectrum this feature peaks 
at considerably shorter wavelengths (« 9.6/im) than ob- 
served for HD179218 or Elias 1 (Hanner et al. 1994). Its 
intensity is small in comparison to the underlying con- 
tinuum, with a ratio of the peak - to - continuum flux 
of about 1.4. The FWHM as derived by the ISOCAM 



CVF images or the ISOPHT multi-aperture sequence of 
HR6000 and HR5999 is typical for point sources. 

3.5. T Cha 

T Chamaeleonis is a T Tauri star of spectral type G8 (Al- 
cala et al. 1993) and the coolest star in our sample. It 
is superposed on a dark cloud DC 300.2—16.8 between 
the star forming Chamaeleon I and II clouds. Although 
a priori the later spectral type suggests that it is physi- 
cally more difficult to excite transiently heated particles, 
we considered this star suitable for our sample as its pre- 
sumed parent cloud DC 300.2—16.8 shows an unusually 
high amount of excess radiation in the mid-infrared (Lau- 
reijs et al. 1989). 

The Hipparcos distance to T Cha is 66j^2 P c - This 
is considerably closer than the estimated distances to the 
Chamaeleon I (160±15) and II (178±18)pc clouds (Whit- 
tet et al. 1997). On the other hand, Terranegra et al. 
(1999) associate T Cha to a moving group of pre-main 
sequence stars which have the same radial velocity as the 
CO gas in DC 300.2-16.8. They note, however, that de- 
spite sharing a common proper motion with the group, 
T Cha has a deviating radial velocity measured by Covino 
et al. (1997). Therefore, unfortunately, it is not possible 
to make any firm statements of the possible physical con- 
nection between T Cha and DC 300.2-16.8. 

The spectrum of T Cha is featureless, with no indica- 
tion of IEB or silicate emission (Fig. |l5|). The continuum 
we measure is about a factor 4-5 larger than the expected 
photospheric flux (about 0.2 Jy at 6 /im), possibly due to 
a circumstellar disk (Alcala et al. 1993). Given the lack of 
dust features it is not a surprising result that there is no 
extended mid-infrared emission around T Cha. Also the 
ISOPHT multi-aperture sequence is for apertures below 
23" consistent with the point source HR7127. It is pos- 
sible that a G8 star cannot excite extended mid-infrared 
emission. However, given the fact that the connection to 
DC300.2-16.8 is uncertain, we cannot exclude that the lack 
of extended IEB emission is simply due to the absence of 
low-density circumstellar matter associated with the star. 

4. Conclusions 

We find that the mid infrared spectrum of the Her- 
big Ae/Be stars: HD 97048, HD 97300, TY CrA and 
HD 176386, is dominated by IEB and show no signature 
of silicate emission. All those stars have an extended halo 
on scales of up to a few thousand AU. Two of the stars 
(TY CrA and HD 176386) have a secondary peak of emis- 
sion at distances of about 1000 AU from the main compo- 
nent. A secondary peak is found also in HD97300 and in 
addition, there is emission in an elliptical ring structure 
at a much larger scale (~ 7500 AU from the star). 

The IEB ratios in the secondary peaks show differences 
when compared to those of the central regions. This can be 
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explained as being due to differences of the excitation of 
transiently heated particles and variations of the column 
density (Siebenmorgen et al. 1998). 

The other targets of our sample are at distances simi- 
lar to the IEB dominated objects. They remain point like 
at the spatial resolution of ISO (~ 5" at 6/zm). Two of 
those stars (HR6000 and T Cha) show featureless mid IR 
spectra while the other objects (HR5999 and HD 179218) 
have a pronounced silicate emission bump. Because sili- 
cate grains are much larger particles than the IEB carri- 
ers (Siebenmorgen et al. 1992), they arc not transiently 
heated and therefore must be located close to the star, 
most likely in a disk component. This interpretation is in 
agreement with the compact emission detected. The mid 
IR spectrum of HD 179218 shows both silicate emission 
and IEB; nevertheless it is not extended. Consequently 
the detection of IEB alone is not a clear indicator to as- 
sume an associated extended emission structure of more 
than a few hundred AU. 

Although the ISOCAM CVF has a known stray-light 
component we could confirm the deduced sizes (FWHM) 
by ISOPHT multi-aperture photometry. We find that the 
total flux is in the ISOSWS spectra of the extended sources 
systematically underestimated. This is expected because 
the ISOSWS flux calibration assumes a point source. For 
point like sources we derive for both ISO instruments a 
set of consistent photometric flux density values. 
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aperture photometry on the CVF images. 
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Fig. 10. ISOPHT multi-aperture sequence at 7.3 on 
HD 176386. For apertures larger than 80" the adjacent 
source TY CrA dominates the measured signals (Fig. 0). 




Fig. 9. Mid infrared spectrum of the second component 
of HD 176386. 



Fig. 11. ISOPHT multi-aperture sequence, normalised 
flux at 7.3 /itm on a point source (HR 7127) and on the ex- 
tended adjacent sources TY CrA and HD 176386. The full 
and dashed line are for TY CrA and HD 176386, respec- 
tively. The FWHM is deduced as in Fig. 4 and determined 
to 13 ± 3" (TY CrA) and 9 ± 3" (HD 176386). 
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Fig. 12. Mid infrared spectrum of HD179218. The CVF Fig. 14. As Fig. |lj for HR5999. 
scan is shown as solid line and the ISOSWS spectra by 
the dotted line. 
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Fig. 13. Mid infrared spectrum of HR6000 together with 
a black-body spectrum of 13000 K (solid line). 
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